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In-Gap Spin Excitations and Finite Triplet Lifetimes in the Dilute Singlet Ground 

State System SrCua^Mg^BOs^ 
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High resolution neutron scattering measurements on a single crystal of SrCu2- a: Mg a: (B03)2 with 
x~0.05 reveal the presence of new spin excitations within the gap of this quasi-two dimensional, 
singlet ground state system. Application of a magnetic field induces Zeeman-split states associ- 
ated with S=l/2 unpaired spins which are antiferromagnetically correlated with the bulk singlet. 
Substantial broadening of both the one and two-triplet excitations in the doped single crystal is ob- 
served, as compared with pure SrCu2(BC>3)2 . Theoretical calculations using a variational algorithm 
and a single quenched magnetic vacancy on an infinite lattice are shown to qualitatively account for 
these effects. 

PACS numbers: 75.25.+Z, 75.40.Gb, 75.40.-s 
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Quasi-two dimensional quantum magnets which dis- 
play collective singlet or spin gap behavior are very top- 
ical due to the novelty of their ground states^ and their 
relation to high temperature superconductivity in the 
copper oxides. There are relatively few such materi- 
als, and crystal growth difficulties have further limited 
their study in single crystal form. SrCu^BOa^ is es- 
tablished as a realization of the two dimensional Shastry 
Sutherland model 2 for interacting S=l/2 dimers 3 ^. It 
is comprised of well separated layers of Cu 2+ , S=l/2 or- 
thogonal dimers arranged on a square lattice. The ma- 
terial crystallizes into the tetragonal space group I42m 
with room temperature lattice parameters of a=8.995 A, 
c=6.649 A£. 

SrCu2(BC>3)2 has been well studied by an array of 
experimental techniques, which show it to possess 
a non-magnetic ground state. In particular earlier 
neutronSi^ and ESR spectroscopyiSiii have estab- 
lished the leading terms in its microscopic Hamiltonian: 
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where J is the exchange interaction within the dimers 
and J' is the exchange interaction between S=l/2 
spins on neighboring dimers. Subleading Dzyaloshinskii- 
Moriya interactions weakly split the three triplet modes 
even in zero applied magnetic fieldZ*&2*i£. 

Both the exchange interactions are antiferromagnetic 
and their ratio x=^j has been estimated between 0.68 
and 0.60, with more recent refinements being smaller. 
Theoretically, such a quantum magnet is known to pos- 
sess a singlet ground state so long as the ratio, x, of inter 
to intra-dimer antiferromagnetic exchange is sufficiently 
smalU 3 .. All of these estimates place SrCu2(B03)2 on 
the low side of the critical value of x at which a quan- 



tum phase transition occurs between a four sublattice 
Neel state and a collective singlet state. 

In finite magnetic field, much interest has focused on 
a finite magnetization which develops at fields beyond ~ 
20 T, wherein the lowest energy of the three triplet states 
has been driven to zero energy 4 * 1 ^ 1 ^ 1 ^. The magnetic 
field acts as a chemical potential for the triplet density 
within the quasi-two-dimensional planes. Magnetization 
plateaus ensue at higher fields, corresponding to Bose 
condensation of the triplets at certain densities. 

While pure SrCu2(B03)2 has been well studied, there 
is little information on this quantum magnet in the 
presence of dopants, and no reports on doped single 
crystals. This problem is very interesting by anal- 
ogy with the remarkable properties of doped quasi- 
two dimensional Mott insulators and high tempera- 
ture superconductivity^. SrCu2 (BC>3)2 is itself a Mott 
insulator and the theory of doped Mott insulators 
on the Shastry-Sutherland lattice shows the possibil- 
ity of several different superconducting phases as a 
function of dopingi2ii2i2S. Several doping studies of 
SrCu2(B03)2 have been reported on polycrystalline sam- 
ples wherein Al, La, Na, and Y substitute at the Sr site 
and Mg substitutes at the Cu site2i. 

In this Letter we report high resolution time-of-ffight 
neutron scattering measurements on large single crystals 
of SrCu 2 - I Mg :c ( 11 B03)2and SrCu 2 ( n B0 3 )2 . These 
measurements show that doping of the magnetic Cu 2+ 
site with non-magnetic, isoelectronic Mg 2+ at the 2.5 % 
level introduces new magnetic excitations into the singlet 
energy gap, and gives a finite lifetime to all three single 
triplet excitations, while also substantially broadening 
the two triplet bound state. 

Two single crystal samples, SrCu 2 ( 11 B03) 2 and 
SrCu2- K Mg a; ( 11 B03) 2 , were grown by floating zone im- 
age furnace techniques at a rate of 0.2 mm/hour in an 2 
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FIG. 1: (a) The two left hand panels show neutron scat- 
tering data from SrCu2(B03)2 at T—2K. The scattering has 
been integrated along L and we show data in a magnetic 
field of zero and 7 T. The right hand panels show the 
same data for SrCu2- a; Mg a; (B03)2 . The energy axis is 
on a logarithmic scale. The Zeeman-split S=l/2 level in 
SrCu 2 _ a; Mg a; (B03)2 in H=7 T at g^ s H=0.8 meV is indicated 
in panel d). 
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FIG. 2: The right hand panels show inelastic neutron scat- 
tering data for both SrCu2(B03)2 and SrCu2-iMgi (BOs)2 in 
H=0 (b) and H=7 T (d). This data has been integrated 
in Q between H=l and H=3 and over all L. The inten- 
sity is plotted on a logarithmic scale. Clearly the one and 
two triplet excitations are considerably broader in energy in 
SrCu2- a; Mg a ;(B03)2 as compared with SrCu2(BOs)2 , and in- 
gap states are introduced on doping. The left hand panels 
show the numerical calculation for S zz (Q,lj) + S yv (Q,u), 
using J=76.8 K and J'/J=0.62, which accounts for many of 
the qualitative features observed in SrCu2- a; Mg a; (B03)2 ■ 



atmosphere. The crystals were of almost identical cylin- 
drical shape, with approximate dimensions of 4.5 cm in 
length by 0.6 cm in diameter. These samples were grown 
using to avoid the high neutron absorption cross 

section of natural boron. The pure SrCu2( 11 B03)2 single 
crystal was the same high quality single crystal studied 
previously 9 . Time-of- flight neutron scattering measure- 
ments were performed using the OSIRIS spectrometer^ 
at the ISIS Pulsed Neutron Source of the Rutherford 
Applcton Laboratory. OSIRIS is an indirect geometry 
time-of-flight spectrometer which employs an array of 
pyrolytic graphite monochromators to energy analyse 
the scattered neutron beam. The data was collected with 
the spectrometer configured to use the 004 analysing 
reflection afforded by pyrolytic graphite. In this con- 
figuration, only those neutrons with a scattered energy 
of 7.375 meV are Bragg reflected towards the detector. 
The [H,0,L] plane of both crystals was coincident with 
the horizontal scattering plane, and the samples were 
mounted in a 7 T vertical, [0,K,0], field magnet cryostat. 

Figure 1 shows representative time-of-flight neutron 
scattering data, taken at T=2 K and H=0 and 7 T for 
SrCu 2 (B0 3 ) 2 (a and c) and for SrCua-zMg^BOsja (b 
and d). This data was integrated along L, in which direc- 
tion the spin excitations show little dispersion 9 .. Note the 
logarithmic energy and intensity scales, chosen to draw 



attention to the details of the in-gap excitations seen in 
SrCu2- a: Mg a ;(B03)2 • The splitting of the triplet exci- 
tations near 3 meV on application of the 7 T magnetic 
field is clear. In finite field, weak dispersion of the triplet 
excitations as a function of wavevector H is seen, and 
this has been attributed to subleading terms in the spin 
Hamiltonian - terms other than those in Eq. 1. There 
are several qualitative features evident on examination of 
this data. The one triplet excitations show significantly 
greater breadth in energy in SrC^-^Mg^BOa^ than in 
SrCu2(B03)2 ■ In addition, application of a H=7 T mag- 
netic field gives rise to an inelastic peak at g/isH =0.8 
meV in H=7 T, which is centered around [H~ 1.4, 0, 0] in 
Q-space, but extends in wavevector H to almost [H^3, 0, 
0] . Both of these features are discussed at length below. 

Figures 2b and 2d show the same experimental data 
as in Fig. 1, now integrated in wavevector along L and 
also in H between H=l rlu and H=3 rlu, and plotted as 
a function of energy. The data in a magnetic field of 
and 7T is shown in Figs 2b and 2d, respectively. The 
extra breadth in both the single triplet excitations and 
the two triplet bound states above it is clear. Broad in- 
elastic peaks appear within the gap, which possess little 
magnetic field dependence, indicating a longitudinal na- 
ture. Sharp, field-induced inelastic scattering at tko ~ 
0.8 meV is also evident. 
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Numerical calculations have also been carried out us- 
ing a new variational approach-- to solve the model 
of a single quenched impurity on the two dimensional 
Shastry-Sutherland lattice. This method generates 
a variational space by successively applying the off- 
diagonal parts of the Hamiltonian, Eq. 1, on the start- 
ing approximation for the single impurity ground state, 
which consists of a free spin 1/2 neighboring the impu- 
rity site, embedded within a dimer background. The 
resulting small spin polaron structure and exponential 
growth of the variational space with each iteration guar- 
antee good convergence of the spin polaron ground state, 
as well as for the lowest energy excited states. Full de- 
tails will be given separately 23 . For energies below ~ 
3 meV the method provides accurate and converged re- 
sults for both the longitudinal, S zz (Q,u>), as well as the 
transverse, S ya (Q,w) components of the dynamical spin 
structure factor. These are compared directly to the 
neutron scattering experiments in Fig. 2. 

Figures 2a and 2c show the calculated S ZZ (Q, u>) + 
S VV (Q, u>), integrated over the same range of wavevectors 
as the experimental data, and at magnetic fields of 
(a) and 7 T (c). The comparison between theory and 
experiment is qualitatively good. The numerical results 
confirm that quenched magnetic vacancies induce in-gap 
states, substantial spectral weight below the zero field 
gap energy of ~ 3 meV. A component of these in-gap 
states show little magnetic field dependence, and appear 
in the S ZZ (Q, u>), longitudinal channel. The calculation 
also shows the g//sH transverse, S vv (Q,uj), excitation 
in finite magnetic field. Furthermore, the calculation^ 3 , 
allows a determination of the spatial distribution of the 
spin polaron S=l/2 degree of freedom and its low lying 
excited states, and these can guide the interpretation of 
the Q-dependence of the field-induced g/^sH transverse 
spin excitation observed in the experiments. 

Figure 3 shows the Q dependence of the scattering 
around g/isH=0.8 meV in an applied magnetic field of 7 
T in both SrCu 2 (B0 3 ) 2 (a) and SrC^-xMg^BC^a 0). 
This scattering is integrated in wavevector over L, and 
in energy between 0.6 meV < tuv < 1 meV and is 
shown in a magnetic field of both and 7T. Figure 3a 
shows the absence of a magnetic field induced signal in 
SrCu2(B03)2 within this energy range. Figure 3b shows 
a clear field induced signal in SrC^-^Mg^BOs^ which 
peaks at [H~1.4, 0, 0], but extends out to almost [H~3.0, 
0, 0]. This field induced scattering is attributed to Zee- 
man split S=l/2 states associated with the S=l/2 mo- 
ment in a dimer whose partner site is occupied by a 
quenched, nonmagnetic Mg 2+ ion. This field induced 
inelastic scattering is very similar to that associated 
with end states in Haldane spin chains, such as occur 
in Y2BaNii_ :E Mg x 05 24 . In this case, quenched, non- 
magnetic Mg 2+ ions produce finite S=l magnetic chains. 
Spin 1/2 degrees of freedom arise at the end of finite 
chains of S=l magnetic moments, as one of the two ef- 
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FIG. 3: Q-scans are shown at T=2 K for both 
SrCu 2 (B0 3 )2 (a) and SrC^-^Mg^BOa^ (b) which inte- 
grate in energy between 0.6 meV and 1 meV. This energy 
range captures the Zeeman energy g/isH=0.8 meV appropri- 
ate to g=2 and H=7 T. (b) shows a comparison between the 
experimental Q-dependence of this scattering with the cal- 
culated form described in the text. The field-induced in-gap 
state in SrCu2-iMg 2! (B03)2 is absent in SrCu2(BC>3)2 . 



fective S=l/2 degrees of freedom making up the S=l 
moments lacks a partner with which to form a singlet. 
Such excitations occur at an energy of g/isH in finite 
field, and display a wavevector dependence which indi- 
cates antiferromagnetic correlations into the collective 
singlet of the chain segment. 

The wavevector dependence of the magnetic field- 
induced spin excitation at 0.8 meV can be attributed 
to the structure of the spin polaron 2 ^, whose ground 
state possesses strong antiferromagnetic correlations 
with neighboring dimers, transverse to the dimer con- 
taining the impurity site. To first non-trivial order in 
J'/ J, we obtain an analytical expression for the square 
of the matrix element for this transition: 

I* <x (a 2 - a\/2) e "'(H±K) + 2 ^7> aifl2 x 

sin(?7(H T- K)) sin(7r(H qp K)) - 2<S £12(23 X 

cos(?7(H =p K)) cos(tt(H t K)) ^ (2) 

where ax, a,2, and CI3 are the weights of the polaron vari- 
ational wavefunction computed for J'/ J = 0.632^. 77 = 
0.72 accounts for microscopic distances in SrCu2(B03)2- 
The ± sign in Eq. distinguishes between the two 
nonequivalent impurity positions within the unit cell. 
Random doping therefore implies I cx I + + I~ , which 
we show, multiplied by the square of the magnetic form 
factor for Cu 2+ , in Fig. 3b along with measurements on 
Q-dependence of this scattering in SrCu 2 - a; Mg 2 ; (603)2- 
While the agreement between the calculation and exper- 
iment is not perfect, the calculation captures the general 
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FIG. 4: (a) Cuts of the data simulating constant-Q scans at 
(1.5,0,L), integrated along L, are shown for SrCu2(BC>3)2 and 
SrCu2- :E Mg a ;(B03)2 • Resolution convoluted fits to the data 
are shown as the solid lines and the description of the data 
is excellent. Such fits allow us to extract the triplet exci- 
tation lifetimes (F), and which are shown in (b), (c), and 
d) for wavevectors (1.5, 0, F), (2, 0, F) and (2.5, 0, F) 
respectively, as a function of temperature. Even at the 
lowest temperatures, we observe finite triplet lifetimes in 
SrCu 2 -*M g;c (B03)2. 



Q-dcpendence of the excitation. 

Figure 4a shows representative data with accompany- 
ing fits used to extract the lifetimes of the one triplet 
excitations as a function of doping and temperature 
in a magnetic field of 7 T. This data is integrated in 
Q along F, and over a narrow range in wavevector H 
around H=1.5, and these cuts approximate constant-Q 
scans. Similar analysis was performed around wavevec- 
tor H=2.0 and 2.5 to look for variation in the one-triplet 
lifetimes as a function of Q along [H, 0, 0]. The data 
was fit to the sum of three identical damped harmonic 
oscillators (DHO): 
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5(Q,w,T) = x(Q,T) 



1 - exp {—fe 



where x(Q, T) is the momentum-dependent suscepti- 
bility, while the second term takes into account de- 
tailed balance. The renormalized DHO frequency, Qq, 
has contributions from the oscillation frequency, ujq.t, 
and the damping coefficient, Tq^, and is given by: 

n Q,T = W Q,T + r Q,T- 

Equation 3 was convoluted with an appropriate reso- 
lution function and fit to the data. The results are shown 
in Figs. 4b), c), and d), for wavevectors H=1.5, 2.0, and 
2.5 respectively. The extracted lifetimes show little or 
no systematic variation with wavevector H, but a finite 
one-triplet lifetime is observed in SrC^-zMg^ (1303)2 at 
all temperatures, in contrast to SrCu2(B03)2 , where the 
low temperature one-triplet lifetimes are very long, com- 
pared with the resolution of the spectrometer. In both 
SrC^-zMg^BOa^and SrCu 2 (B0 3 )2 , the thermal de- 
struction of the collective singlet ground state near ~ 10 
K is characterized by a rapid decrease in the one triplet 
lifetime (1/r) on warming, with little or no softening of 
the one triplet excitation energies. 

As mentioned previously, the theoretical results for 
S zz (Q, w) + S vv (Q, uj) are not well converged for energies 
of ~ 3 meV and greater. Nonetheless, the additional 
broad spectral weight around the calculated one and two 
triplet energies seen in Fig 2a) and c) is consistent with 
finite triplet lifetimes in the presence of a single quenched 
magnetic vacancy. 

To conclude, new inelastic neutron scattering mea- 
surements on SrCu2- a; Mg a ;(B03)2 with x~ 0.05 show rel- 
atively broad and field-independent in-gap spin excita- 
tions as well as a magnetic field-induced excitation iden- 
tified as a Zeeman-split, spin polaron state. The non- 
magnetic quenched vacancies also give rise to finite and 
measurable lifetimes in the one and two triplet excita- 
tions. 
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